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Abstract. In this paper, we study the bilinear form and the general V-soliton solution for 
a two-component Hunter-Saxton (2-HS) equation, which is the short wave limit of a two- 
component Camassa-Holm equation. By dehning a hodograph transformation based on a 
conservation law and appropriate dependent variable transformations, we propose a set of 
bilinear equations which yields the 2-HS equation. Furthermore, we construct the V-soliton 
solution to the 2-HS equation based on the tau functions of an extended two-dimensional 
Toda-lattice hierarchy through reductions. One- and two-soliton solutions are calculated and 
analyzed. 
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1. Introduction 

The Hunter-Saxton (HS) equation 


Htxx 2 ,KUx 4“ 2,UxUxx T UHxxx — O 5 


( 1 . 1 ) 


was derived as a model for propagation of orientation waves in a massive nematic liquid 
crystal director field |[T1 . The Lax pair, the bi-Hamiltonian structure and local and global weak 
solutions were discussed by Hunter and Zheng El ■ It belongs to the member of the hierarchy 
of the Harry-Dym equation lO, and in a series of papers by Alber and his collaborators, the 
link between billiard solutions and soliton solutions of HS equation were made clear [fflEllbl. 
The HS equation can be regarded as a short wave limit of the well-known Camassa-Holm 
equation lUl [H 


rrit + Unix + 2mUx = 0 , m = K + u — Uxx ■ 


( 1 . 2 ) 


In the present paper, we are concerned with the two-component Hunter-Saxton (2-HS) 
equation 


Htxx 2 kUx “h 2UxUxx T UUxxx — 


(1.3) 
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P? + (pM)x = 0, (1.4) 

where k and a are positive eonstants. The 2-HS equation has attraeted mueh attention in the 
past and it has been studied extensively by many authors ^ [TOl [HI [l3l HH [T5l HH [TTll . 
In a series papers by Wunseh, he studies the loeal and global weak solutions for the periodie 
2-HS equation. The single solitary wave solution was studied in [fT^fTSlfTTll . The loeal well- 
posedness and wave-breaking was studied by Moon and Liu lfT4l . The geometrie property of 
the 2-HS equation was investigated in ifT^ . 

The two-eomponent Hunter Saxton system is a partieular ease of the Gurevieh Zybin 
system pertaining to nonlinear one-dimensional dynamies of dark matter as well as nonlinear 
ion-aeoustie waves (ef. llTSl and the referenees therein). It was also derived as the N = 2 
supersymmetrie extension of the CH equation [fTOl . It is sometimes ealled the generalized 
two-eomponent Hunter-Saxton system if K 7 ^ 0. We simply eall it the two-eomponent Hunter- 
Saxton equation for either K = 0 or K 7 ^ 0 hereafter. The 2-HS equation ean be viewed as 
the short wave limit of the two-eomponent Camassa-Holm (2-CH) equation originating in 
the Green-Naghdi equations whieh approximate the governing equations for water waves 


Ellaoiizii 

mt-\-umx-\-2mUx — Oppx = 0, (1.5) 

P? + (pM)x = 0, (1.6) 

m = K-^u-Uxx, (1.7) 

by using the sealing (t,jc) —)■ (£t,£^^x). To be speeifie, let 

x = et, ^ = (1.8) 

and expand u and p in power series as 

M = £^(mo + £mi-T ■ ■ ■), p = £(po + £pi + ■ ■ ■). (1-9) 

Then we have dt = edi, dx — Substituting these relations together with (11.91) into Eqs. 

(11.5I) - (I1.7I) . we obtain the following PDE for uq and po at the lowest order in £. 

2kmo,^ - - 2.uq iUq '^^ - uqUq^^^^ -b apopo,^ = 0, (1.10) 

Po,T +(poMo)^ = 0. (1.11) 


Writing baek above equations in terms of the original variables t, x, we arrive at the two- 
eomponent Hunter-Saxton equation (ll.3l) - (ll.4l) . 

Both the two-eomponent Camassa-Holm equation and two-eomponent Hunter-Saxton 
equation are integrable in the sense that they admit the Eax pair 

( 1 . 12 ) 

% = -(^-^ + u^^'x+\ux^. (1.13) 

where m = k-\- u — Uxx for the 2-CH equation (ll.5l) - (ll.7l) and m = K — Uxx for the 2-HS 
equation (11.3I) - (I1.4I) . The eompatibility eondition gives the 2-CH and 2-HS equations, 
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respectively. Therefore, they both have infinite numbers of conservation laws and bi- 
Hamiltonian structures. 

When K 0, the 2-HS equation (ll.3l) - (ll.4l) is invariant under the following scaling 
transformations 

m^Em, p-)-£p, K-)■ £K, a£^a, (1-14) 

Therefore, without loss of generosity, we can fix one of the values for either K or o. 

In spite of many studies regarding the integrable 2-HS equation, as far as we are aware, 
its the two- and general A^-soliton solution are not known yet. Therefore, the goal of the 
present paper is to find multi-soliton solution to 2-HS equation by using Hirota’s bilinear 
method [l22l . The remainder of the present paper is organized as follows. In Section 2, 
through the hodograph transformation implied by the second equation in 2-HS equation, we 
propose a set of bilinear equations for 2-HS equation, and present its one- and two-soliton 
solutions by Hirota’s method. In Section 3, Starting from bilinear equations of an extended 
two-dimensional Toda-lattice hierarchy, together with their Casorati determinant solution, we 
derive the general multi-soliton solution through a period 2-reduction. A detailed analysis 
for one- and two-soliton solutions shows that there are either smooth or loop soliton solution 
depending on the choices of parameters. Section 4 is devoted to some concluding remarks. 


2. Hodograph transformation and bilinear equations for two-component 
Hunter-Saxton equation 


The second equation (11.41) of 2-HS equation represents a conservation law, from which we 
can define a hodograph transformation, or sometimes called a reciprocal transformation 

dy = pdx — pudt ^ ds = dt. (2.1) 

It then follows that 


dx — pdy, dt — ds pudy , 
or 

dy — p dx , d^ — df -j- udx . 
Using (12.21) . Eq. (11.41) can be rewritten as 
Pi 3“ P Uy = 0 1 
or, equivalently. 


= u 


y 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 


In view of the dependent variable transformations for the Camassa-Holm equation and short 
wave limit of the Camassa-Holm equation lf23l[24ll . the following transformation 

,2 


w = -(lng)ii, p = 


8 ^ 

fh 


( 2 . 6 ) 
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is suggested, where /, g and h are usually tau funetions. Substituting the transformation (12.61) 
into (12.51) and integrating with respeet to s, we arrive at the following bilinear equation 

\DyDs-l^g-g^-fh, (2.7) 

by taking an appropriate integration eonstant. Here D is ealled Hirota /^-operator defined by 

d d / d d 


D^:D;y-g = 


ds ds' 


f{y,s)g{y ,s)\y=y^s=s'- 


dy dy' 

We remark here that the bilinear equation (12.71) ean be viewed as the one for two-dimensional 
Toda-lattiee, whieh will be made more elear in the subsequent seetion. Note that the first 
equation of 2HS takes an alternative form 

mt + umx + 2mux-0ppx = 0, ( 2 . 8 ) 


m = 2 — u 


XX 


(2.9) 


where k is set to be 2 due to the invariant property (11.141) of the 2-HS equation (ll.3l) - (ll.4l) . 
An substitution of (12.31) into (12.81) and (12.91) yields 

ms-\-2mpuy — (5p Py = 0, 


m = 2 + p(lnp)y^, 

respeetively. Eqs. (I2.10I) - (I2.11I) imply the following two bilinear equations 


DyDs -|- + sf^Dy ) / ■ h — 0,, 


D, + l\f-h = g^ 


1 

\/a' 

In summary, we have the following theorem. 
Theorem 2.1. A set of bilinear equations 

^-DyD,-]^g-g^-fh, 

1 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


-^D, + ljf-h = g\ 

DyDs + ~^J^Ds + s/oDy ) / ■ = 0 , 


(2.14) 

(2.15) 

(2.16) 


give the two-component Hunter-Saxton equation ( li.3D -( l7~?l) with K = 2 through the dependent 
transformation 

A 


“ = -(lng)ii, p = 


g 


fh' 


(2.17) 


and the hodograph transformation 

x = y-(lng)^, t = s. 


(2.18) 
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Proof. Dividing both sides of Eq. (l2.14l) by g^, we have 
n ^ 1 

{\ng)ys-l = -^. 

Dividing both sides of Eqs. (12.151) and (12.161) by fh, we arrive at 




and 


+ 1 = ^, 
. fh' 


f 




(ln//.V+ In^ + _2 = 0, 




/ 


\/o 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


respeetively. Multiplying Eq. (l2.19l) by 2, then subtraeting it from Eq. (12.211) . one yields 


fh 


6 / ys 


f 




+ In^ +^/^ int + ^ = 2 ^ ■ 




/ 


Jh 


s/o 




( 2 . 22 ) 


By referring to Eq. (12.201) . it ean be easily eheeked that 

/IN 1 /I N -1 

^ = -(lng),, = M, — = 1 - (Ing)y, = p 

whieh realizes the hodograph transformation (12.11) defined previously. Note that Eq. (12.191) 
ean be written as 

, 

or equivalently 

(Inp)i — pUy = Ux ■ 

Combining Eq. (12.221) with Eq. (12.201) , we have 


(2.23) 


(2.24) 


+ 


r 


(2.25) 


(2.26) 


2 + p(lnp)y, = \/op2 I (In 

Eet us define 

m = 2 + p(lnp)y^, 
then it immediately follows 

m = 2-Uxx, (2.27) 

by using Eq. (12.241) and the eonversion relation (12.31) . Differentiating Eq. (12.251) with respeet 
to we obtain 


m, 


- 2^/opp, (In 


+ 


y 


^/o 


+ ^/op^ (In 


/ 


ys 


2 + p(lnp) / f 

= 2p,- + V^p2 In L 

P \ h 

= 2m{\np)s + op^Py 


ys 


= -2mux + (5ppx. 


( 2 . 28 ) 
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Here equations (12.201) . (12.241) and (12.251) are used. 

Finally the hodograph transformation (12.31) eonverts Eqs. (12.241) and (12.281) into 

{dt + Mdx)p = -pMx, (2.29) 

and 


{dt-\-udx)m~—lmUx-\-<5pPx^ (2.30) 

respeetively. The above two equations are nothing but the 2-HS equation with m = 2 — Uxx- □ 


Remark 2.2. Substituting (12.151) into (12.151) . one obtains 

[DyDs + s/oDy -2)f-h = 2g^. (2.31) 

If we assume / = /z as a —)■ 0, then Eq. (12.311) eonverges to 

(^^DyDs-\^f-f = 2g\ (2.32) 

It is worth pointing out that Eqs. (12.161) and (12.311) are aetually the bilinear equations for 
period 2-reduetion of two-dimensional Toda-lattiee, whieh yield the sine-Gordon equation 
and further the short-pulse equation and the short wave limit of the CH equation. 


In the last, we proeeed to finding the one- and two-soliton solutions for 2-HS equation 
(ll.3l) - (ll.4l) by Hirota’s perturbation method ll22l . 

One-soliton solution: To find one-soliton solution, we assume 

g=l+e*?'+“^, /=l+ce^+“^ /^=l+Je^^+“^ (2.33) 

and substitute them into the bilinear equations (I2.16I) - (I2.15I) . As a result, we obtain the 
following algebraie relations 

mk = 2 —{c-\-d), cJ = 1, 

—;=(c — d) =2 — {c-\-d), 

Va 

(Ok -T \/ok{c — d) = 2(c — J) — 4. 


By ehoosing c as a free parameter, we yield the one-soliton solution determined by 


k = 


1 - 1 

(0 = vo- 


d = 


1 


"v/o c ’ '' 1 -|-c c 

The tau funetions found here leads to the one-soliton solution of the form 


o 

~4 


1 —c 
l-fc 


seeh ( -{y-vs) ) , 


p=|l+(i^,ech^ 

' 4c 


;(y-v5) 


(2.34) 

(2.35) 


■ = y -h \/0' 


1 


1 —c 


l+c V 


- 1 


t = s. 


(2.36) 


The amplitude of the soliton is a(l — c)^/(4(l -l-c)^), and the veloeity in the (y,^) plane is 
V = oc/{I -\-cY. Erom Eqs. (I2.34l) - l234l) . it is obvious that the eonditions of c > 0 and c 7 ^ 1 
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must be imposed in order to assure a non-singular soliton solution. It is interesting to note that 
there is an upper bound for the amplitude of solitons, which is o/4. When c —)■ 0 or c —)■ 0 °, 
the amplitude of the soliton approaches this upper bound, meanwhile, the velocity approaches 
zero. When c ^ 1, the amplitude approaches zero, whereas, the velocity approaches to the 
maximum of o/4. An example with c = 7 and o = 4 is illustrated in Fig. 1. 




Figure 1. An example of soliton solution for 2-HS equation: (a) x-u profile; (b) x-p profile. 


Two-soliton solution: Based on the obtained one-soliton solution, we assume the tau 
functions of two-soliton solution are of the form 


g = I gkiy+c)is ^k2y+(i)2S ^^^g{ki+k2)y+{(^i+C32)s ^ 


(2.37) 

(2.38) 

(2.39) 


where 




(0 = 


1 




1 , 2 ). 


y/oCi 1 +C, 

After some tedious calculation, we obtain two-soliton solution determined by 


an = 


Cl -C2 Y 

1-C 1 C 2 / 


ci2 = ciC2ai2, dn = c^^C2^an- 


We illustrate two-soliton solution interaction in Figs. 2-4. As shown in Fig.2, initially, we 
choose one soliton with c = 2.0 located on the left, and one soliton with c = 7.0 located on 
the right. They undertake a collision later on. The collision is elastic without any change 
of shape but a phase shift after the collision which is shown in Fig. 2-4 for f = 0,30,60, 
respectively. 
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Figure 2. Two-soliton interaction at f = 0; x-u profile; (b) x-p profile. 




Figure 3. Two-soliton interaction at f = 30: x-u profile; (b) x-p profile. 




Figure 4. Two-soliton interaction at f = 60: x-u profile; (b) x-p profile. 
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We start with a Casorati determinant solution for the extended (deformed) two-dimensional 


Toda hierarchy 


where 


— 



(n) 

¥i 

(«) 

¥2 

vr" .. 
>,,<"+» .. 

(n+N^l) 

■ Y1 

(«+yv-i) 

Y2 

5 

( 3 . 1 ) 

In) 

Tn 

(17+1) 

Vn 

Y^V 



= ai,\ 

(PT'-bTe'^ 

+ai, 2 {q 7 ^ - 

bYe^‘ , 

( 3 . 2 ) 


with 




-x\ -\- 


Pi 


(Pi -b) 


-X2 Pi Pi ^X-2 + , 


(3.3) 


1 


qi -b 


x\ -\- ■ 


1 


;X2 + qi x-i+qi X-2+r\i0. 


(3.4) 


{‘li ^-bY' 

The following lemma gives bilinear equations satisfied by the above tau-functions presented. 
Lemma 3.1. The tau functions x„ given by t l3. 1 D -t lX?!) satisfy the bilinear equations 

^ - ■ ■ (3.5) 

(3.6) 

(3.7) 

(3.8) 


Dx\Dx \ ^ J l^n +1 5 

P^xi {P^x-2 )Xn ■ Xfj ■ X^j—j, 

Dx^P^X I'^n ■ ~ Dxi'^n-Vl ' ; 

1 


7^X2(7^x_2 2Z7Z)x_j) 2 I x^ ■ Xfi — {Dx\T^x I 2)Xn^i 'X^—i. 


Proof. The proof will be given by the technique developed by Hirota and Ohta [l22l l25ll . 
Firstly, we can easily check that the following relations 


-v (n) (77— 1 ) 

( 3 . 9 ) 

-V (77) (77 — 2) 

5 x 2 ¥) - ^ 

( 3 . 10 ) 

-v (77) (77+1) , (77) 

( 3 . 11 ) 

-V (77) (77+2) , n.1 (77+1) , ,2 (77) 

ax_2¥/ =Ti + + b^Ti ■ 

( 3 . 12 ) 


For simplicity, we introduce a convenient notation. 



¥ 2 ”'^ ■ 


¥("’ 

■ 



(3.13) 
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by which is rewritten as 

X„ = |n,n + 1, • • • ,n + A^—1|. (3.14) 

Based on above relations (I3.9I) - (I3.12I) . we have the following relations regarding the 
derivatives of x„ 

dxi^n = \n-l,n + l,---,n+N-l\, (3.15) 

(dx_i -Nb)Xn = \n,n+l,--- ,n + N -2,n + N\, (3.16) 

{dx,{dx_,-Nb)-l)Xn=\n-\,n + \,---,n+N-l,n+N\. (3.17) 

(dx _2 — ^b^x—i 3 “bib ')Xfi = |n,n-|-l,‘'',n-l-A^ — 2,n-l-A^-|-l| 

— |n,n+l,---,R + 7/ — 3,n+77—l,n + A^|, (3.18) 

^xi(5 x_2 “ +Nb^)^n = \n- l,n+I,-■ ■ ,n+N-2,n + N-\-l\ 

— \n — l,n-{-l,---,n-\-N — 3,n-\-N—l,n-\-N\. (3.19) 

5x2(^x_i -Nb)Xn = \n-2,n + l,---,n + N-2,n->rN\ 

— |n — l,n,n + 2,---,n+A^ — 2 ,r + 77|, (3.20) 


dx 2 {^x -2 ~ +Nb^)Xn = \n-2,n+l,---,n+N-2,n + N+l\ 

— |n — l,n,n + 2, •••,n+A^ — 2,n + A^+l| 

— \n — 2,n-\-l,---,n-\-N — 3,n-\-N—l,n-\-N\ 

-{-\n — l,n,n-\-2,---,n-\-N — 3,n-\-N—l,n-\-N\. (3.21) 

Then the Pliieker identity for determinants 

|n—l,n+l,---,n + A^ — 2,n+A^| x |n,n+l,---,n+77 — 2,n + A^—1| 

— \n,n-\-l,---,n-\-N — 2,n-\-N\ x |r— l,n+l,---,n+77 — 2,n + A^—1| 

+ |n+l,---,n + A^ — 2 ,r+A^— l,n + 77| x |r — l,n,n + l,---,n + A^ — 2| =0, (3.22) 

gives 

((^xi (^x_i Nb^ l)'Tn X (^x_i Nb^Xn X d^jX^j + 1 = 0 , 


i.e., 

(djci iX/i Xfi'j X/i dx[ Xndx_ j Xfi -\- Xfi -\-1 Xfi— 1 0, 

which is exactly the first bilinear equation (13.51) . Next, let us prove the second bilinear 
equation. A subtraction of the following two Pliieker determinant identities 

|n—l,n+l,---,n + A — 3,n + A—l,n+A| x |n,n+l,---,n+A — 2,n + A—1| 

— \n,n-\-\,---,n-\-N — 3,n-\-N—\,n-\-N\ x |n — 1 ,r + 1,---,r + A — 2,n+A—1| 

+ |n — l,n,n+l,---,n+A — 3,n + A—1| x |n + l,n + 2,---,n + A—l,n+A| =0, 


(3.23) 
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|« — l,n + l,---,n + A^ — 2,n+A^+l| x |n,n + l,---,n + A^ — 2,n+A^—1| 

— \n,n-\-\,---,n-\-N — 2,n-\-N-\-\ \ x \n — l,n-\-\,---,n-\-N — 2,n-\-N—\\ 

+ |R + l,---,n + A^—l,n + A^+l| X |n — l,n,n+l,---,n+A^ — 3 ,r+A^ — 2| = 0, 


leads to 


i.e., 


^xi {^x ^2 + Nb ) x„ X x,i (^x _2 + Nb x 

3“ (^x_i blb'jXfi+l ^ Xn— 1 \ (^x_i blb'jXfi—l 0, 


(3.24) 


2Z?((dx[3x_iXn)x^ 5x_iX/7dx[Xn) 3“ (3x_[Xn-)-i)Xn_i X^-fi (dx_[Xn.-l)—0, 

which is nothing but the second equation (13.61) . On the other hand, the difference of the 
following two Pliicker determinant identities 

|n — 2,n+l,n + 2,---,n + AC — 2,r+AC| x |n,n+l,---,n+AC — 2,n + AC—1| 

— |n — 2,n + l,n + 2, •••,n + AC — 2,n+AC—1| x |n,n+l,n + 2,---,n+AC — 2,n + AC| 

+ |n — 2,n,n+l,n + 2,---,n+AC — 2| x |n + l,n + 2, •••,n+AC — 2,n+AC—l,n + AC| =0, 

(3.25) 

|n—l,n,R + 2,---,R+AC — 2,n + AC| x |n,n + l,---,n + AC — 2 ,r+AC—1| 

— |n — l,n,n + 2,---,n+AC — 2,n + AC—1| x |n,n+l,n + 2,•••,n + AC — 2,n+AC| 

+ |n — l,n,n+l,n +2,•••,«+// — 2| x |n,n + 2,•••,n4-AC — 2,n + AC—l,n+AC| = 0, 

(3.26) 

implies 

^X2 (^x_i Nb')Xfi Xfi ^X2^n ^ (^x_i Nb^Xfi-\- (^xjXn— 1 )Xn4-l Xn—l (^xiX«-|-l) = 0, 


i.e., 

(dx2^x_iX«)Xn ^x2^«^x_iXn (^xiXn+l)Xn—1 3“ X^+l (dxiXn—l) =0. 

To prove the forth bilinear equation, we need the following four determinant identities 
|n — 2,n3-l,---,n3-A(' — 2 ,r3-AC3-1| x |R,n3-l,---,n3-AC — 2,n3-A('—1| 

— |n,n3-l,---,n3-AC — 2,n3-AC3-l| x |n — 2,n3-l,---,n3-AC — 2,n3-AC—1| 

3-|r — 2,R,R3-l,---,n3-AC — 2| x |n3-l,---,n3-AC — 2,n3-A('—1 ,r3-A('3-1| =0, 

(3.27) 

|n — 2,n3-l,---,n3-AC — 3,n3-AC—l,n3-AC| x |n,n3-l,---,n3-AC—1| 

— |n,n3-l,---,n3-AC — 3,n3-AC—1 ,r3-AC| x |n — 2 ,r3-1,---,r3-AC—1| 

3-|n — 2,n,n3-l,---,n3-A(' — 3,n3-A('—1| x |n3-l,---,n3-A('—I,n3-A('| =0, 


(3.28) 
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jn — l,n,n-i-2,---,n-i-N — 2,n-hN-l-lj x |n,n + l,---,n + A^—1| 

— jn,n-hl,---,n-i-N — 2,n-i-N-hlj X |n—l,n,n + 2, •••,n + A^—1| 

+ |n—+ + — 2| X |n,n + 2,---,R+A^—l,n + A^+l| =0, 

(3.29) 

jn — l,n,n-i-2,---,n-i-N — 3,n-i-N—l,n-i-Nj X \n,n-\-l, - ■ ■ ,n-\-N —\\ 

— \n,n-\-\,---,n-\-N — 3,n-\-N—\,n-\-N\ x |«—l,n,n + 2,---,n + A^—1| 
-\-\n—l,n,n-\-\,---,n-\-N — 3,n-\-N—l \ x |n,n + 2,---,R + A^—l,n+A^| =0. 

(3.30) 

Taking an appropriate linear eombination of these four bilinear identities, and substituting into 
the differential relations for tau functions (I3.15I) - (I3.21I) . we arrive at 

(^^2 (^x_2 -\-Nb ) 2)XnXXn {dx_2 2,bdx_^ Nb )Xn X ^X2^n T X (^x_i 

Nb) l)Xn—1 X X,;-|_i Xn—1 (d^i A^(^) l)X;j-|_i + (^x_i Nb)Xn—l X d^jX^-i-i = 0 , 

i.e., 

(^X2^x_2^”)^« 2Z7( Xn)x,2 ^x_i^«^X2^n) 2XnX,j 

)’^n—1 3“ ^x_i'^n+l^xi 1 T ^Xi'^n+l^x_i '^n—1 '^n+l^Xi^x_i 1 T 2Xn-|-lXn—1 = 0, 

which is exactly the fourth equation (13.81) . The proof is complete. □ 

Next, we show the reduction process in order to construct multi-soliton solution to 2-HS 
equation. We start with a period 2 reduction by requiring 

Gi = -Pi, 


which leads to dx_ 2 '^n = cx„. Thus, the bilinear equations (13.61) and (13.81) are simplified into 


bDxi Dx 1 ’^(1 ■ — b)x ] X, 2 ^ 1 ■ X,j_ 1 , 

(3.31) 

and 

{bDx2Dx 1 T 2) X 22 ■ x^; = {Dx^Dx 1 2)x,2-|-i ■ x^— 1 , 

(3.32) 

respectivelv. A substitution of Eq. 13.51) into Eq. 13.311) leads to 

(^x_i T 2.b^ ■ Xfi—i = 2bXfi ■ Xn , 

(3.33) 

while a substitution of Eq. (13.71) into Eq. (13.321) vields 

{J-2x\Dx 1 T 2,bDxi 2)X2;-|-1 ■ Xn—l — 2Xn ■ Xfi. 

(3.34) 

Finallv. substituting Eq. (13.331) into Eq. (13.341). one obtains 

{DxiDx 1 -\-b ^Dx 1 -\-2bDxi)x„-^-i -Xn-i = 0 . 

(3.35) 


Finally, if we assume Xq = g, Xi = / and x_i = h and let xi — y, x_i = 5 and 2b = y/o, Eqs. 
(13.51) . (13.331) and (13.351) are exactly a set of bilinear equations (I2.14I) - (I2.15I) . which derive 2- 
HS equation. In summary, we can provide the Wsoliton solution to 2-HS equation by the 
following theorem 
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Theorem 3.2. The 2-HS equation ^i.3D - di.4D admits the N-soliton solution of the parametric 
form 


g 

u^-(lng)ss, P = 
x = y-(lng)s, t^s, 


(3.36) 

(3.37) 


where g 


To, / = Xi, /z = x_i with the tau functions x„ determined by 


'^n 






(n+yV-l) 

Y 1 

„>+yv-i) 

Y 2 


where 




in+N-\) 

YA? 





n 

e^‘+ai^2 




with 


(3.38) 


(3.39) 


11 / 


Pi - \/^/2 


y-\-Pi ^5 + ^, 0 , 


Pi ^ + \/o /2 


y-Pi ^+11/0- 


Here a/, 2 . ^/o cind rj/o are arbitrary constants. 


(3.40) 

(3.41) 


The proof is obvious based on the previous diseussions. 

To compare with the one- and two-soliton solution, we list the tau functions for one- and 
two-soliton solutions as follows. 

One-soliton: Let ai 1 = ai 2 = 1, based on the above theorem, we have 


2p. 


g = oc 1 -h 




y+2p^ L+^o 


0^/4 


(3.42) 


/ = {P\^- -(^Pi^+ V^/ 2 ^ e^^ , 

PT^ + ^/^/2 


(3.43) 


h= (t*! ^ - Vo/ 2 'j - (Pi^ + Vo/2^ , 

^ pY^ + V^/2 -^y+2pf^+^o 
P \'-^/2 


(3.44) 
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Two-soliton solution: Based on the general A^-soliton solution, the tau funetions for 
two-soliton solution ean be expanded as 


g = 


^ ~ ~ (P 2 ^ + \/^/2 j 

= (Pl' 

+ (pF' + P2')>t‘^*'^' - (pF*-pf‘)«^‘^^' 

^ I _ pr'+pF/ ^i-m+pd±p£!efc-% 


112 


Pi -P2 


Pi -P2 


By defining 


and 


_ Pi +P2 ^ 4 l- 11 l = g^i-m+^io Pi +P2 ^2-r|2 ^ g^2-r|2+^20 

_1 _1 5 _1 _1 

Pi -P2 Pi -P2 


C; = 


Pi 


Vo /2 


ki = 


2 pr' 


p/ + \/o/ 2 ’ p,. ^-a2/4’ 


CO; = 2 p 


I ’ 


(3.45) 


(3.46) 


(3.47) 


we finally ean show that g is exaetly the expression (12.371) obtained in previous seetion. 
Similarly, we have 

- -(p7^ + (p7^ - f p7^ + 


/ = 


Pi 

P2 


2 ; 


2 J. 


2 A 


_^^e^2_fp2-i + ^^gil2 e^^+(p7^ + 4) 


2 A 


P2 


2 


Pi 


2 7 


j _ Pi ^ +P2 ^ Pi ^ ~ V ^/2 Pi ^ +P2 ^ P2 ^ ~ '^ 1 '^ Xr- 

pr ^- P 2 ^ pr ^+\/^/2 pr ^- P 2 ^ P 2^+\/^/2 

fpF^ - \/o/2') (p2 ^ - A^/2^ 


1)2 


,^i+^ 2 -r|i-il 2 


Pi ^ [Pi ^ + \/^/2 


(3.48) 


and 


/i = 


P/ - 


P2 


2 7 


e- 

^ _ (ry-^ + ^^2 


Pr^ + #) 


P2 + 


2 7 


_|_g 1 l 2 


^ _ P/ + P 2 ^ Pi ^ + Y^jl Pi ^ +P 2 ^ P 2 ^ + A^/2 „^7- 

pr^-P2^pr^-v^/2 P\^ -Pi" Pi^ -sf^Y 

fpi ^ + \/^/2) fp2 ^ + \/^/2 ^ 


1 l 2 


,^ 1 +^ 2 -T| 1-112 


(3.49) 


/Pi ^ - \/o/ 2 ) [P 2 ^ - \/a/ 2 ^ 
whieh eompletely reeover the two-soliton solution (I2.38I) - (I2.39I) through (I3.46I) - (I3.47I) . 
Before we end this seetion, three remarks are given. 

Remark 3.3. If we redefine the parameters as follows 

p/^-\/a/2 , 2pri 


Ci = 


Pi ^ + Y^Y 


ki = 


Pi 


■a2/4’ 


( 0 i = 2 p 


I ’ 
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then we recover exactly the one- and two-soliton solutions for 2-HS equation found in 
previous section. 

Remark 3.4. As a —0, the tau function g goes to the one representing A-soliton solution 
of the Hunter-Saxton equation Il23l . This finding is consistent with the fact that the bilinear 
equations for the 2-HS equation converge to the ones for the Hunter-Saxton equation. In this 
limiting case of one-soliton, ci = C 2 = —1.0, thus, / = /? = ! — g = 1 - 1 - 

Under this case, p becomes singular at the peak point and the soliton is actually the cuspon 
type as mentioned in [l23l . 

Remark 3.5. We should point out here that both the 2-CH equation and the 2-HS equation 
become the CH equation and the HS equation, respectively, in the limiting case of a —)■ 0 not 
p ^ 0. It is interesting that the quantity p still exists in the CH equation and the HS equation, 
it is merely decoupled with the dependent variable u. 


4. Concluding remarks 

In the present paper, we proposed a set of bilinear equations for the two-component 
Hunter-Saxton equation through defining appropriate dependent variable transformations and 
hodograph transformation. Based on this set of bilinear equations, we construct one-, two- 
soliton solution to the 2-HS equation. It is interesting that the one-soliton solution to the 2-HS 
equation is either smooth or loop soliton, whereas, the HS equation only admits the cuspon 
soliton solution and the 2-CH equation admits either smooth or cuspon solution. By a period 
2 reduction of an extended two-dimensional Toda-lattice hierarchy, we provide and prove the 
A-soliton solution in Casorati determinant form. The expansion for one- and two- soliton 
solution agrees with the results obtained via Hirota’s perturbation method. 

In a series of work by one of the authors, the integrable discretizations for a class of 
integrable PDEs with hodograph transformation such as the short pulse equation, the CH 
and its short wave limit, i.e., the HS equation, the reduced Ostrovsky equation and multi- 
component short pulse equation have been developed [[26l [24l |23l |271 |28l, and some of 
them have been successfully used as an integrable self-adaptive moving mesh method for 
the numerical simulations of these PDEs [|2^ |30l . It is a further topic for us to explore the 
integrable discreitzation of the 2-HS equation. 
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